Resonance Patterns of an Antidot Cluster: Prom Classical to Quantum Ballistics 
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We explain the experimentally observed Aharonov-Bohm (AB) resonance patterns of an antidot 
cluster by means of quantum and classical simulations and Feynman path integral theory. We 
demonstrate that the observed behavior of the AB period signals the crossover from a low B regime 
which can be understood in terms of electrons following classical orbits to an inherently quantum 
high B regime where this classical picture and semiclassical theories based on it do not apply. 
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I. INTRODUCTION 

In recent years the physics of electrons in semiconduc- 
tor systems of reduced dimensionality has been studied 
extensively at high magnetic fields where the quantum 
mechanics of Landau levels gives rise to the quantum 
Hall effect, and also at very low magnetic fields where 
many phenomena can be understood in terms of electrons 
following classical ballistic trajectories pi. The crossover 
between these two very different regimes is poorly under- 
stood; identifying and performing relevant experiments 
has been difficult and simple theoretical models that cap- 
ture the essential physics have not been formulated. 

A series of experiments carried out recently in the 
crossover regime have yielded intriguing results . The 
device studied (see the left inset of Fig.[l]) was a quan- 
tum wire defined in a two-dimensional electron gas and 
containing a pair of antidots, i.e., regions from which 
electrons are excluded, in close proximity to each other. 
As a function of magnetic field B, the conductance G 
of this device showed a dip associated with the trap- 
ping of electrons in a closed orbit encircling the antidot 
pair H . Superimposed on this dip were Aharonov-Bohm 
(AB) conductance oscillations due to quantum interfer- 
ence associated with the multiply connected device ge- 
ometry. Plot (a) at the bottom of Fig.|l| is an example 
of such an experimental conductance trace. Within ex- 
perimental uncertainty, the period of the AB oscillations 
was independent of B, a phenomenon that has been 
difficult to understand since the area enclosed by a clas- 
sical cyclotron orbit decreases as B increases. One might 
expect this to result in an AB period that increases with 
B. [Q Another intriguing experimentally observed phe- 
nomenon that will be reported here is that when the ge- 
ometry of the antidots is varied smoothly by varying the 
voltage on the gates that define them, an additional AB 
conductance peak can, under certain conditions, appear 



abruptly, as if the electronic orbit suddenly enclosed an 
extra quantum of magnetic flux. As the gate voltage 
is varied further, the sequence of AB resonances (as a 
function of B) soon reforms and becomes regular again. 
Such dislocations in the pattern of AB resonances are 
observed when more than one transverse mode is trans- 
mitted through the constrictions between the antidots 
and the edges of the quantum wire. The independence of 
the AB period of the value of B is observed both in this 
regime and when only a single mode is transmitted. 

This is clearly an ideal system for comparing classical 
and quantum behavior. In this article we account for 
its surprising behavior with the help of computer simu- 
lations and Feynman path integral theory. Our quantum 
simulations reproduce the experimentally observed con- 
ductance dip and AB resonances. We show the latter to 
be due to resonant transmission (as distinct from reso- 
nant reflection) of electrons. We demonstrate that the 
experimentally observed constancy of the AB period is 
an innately quantum phenomenon that signals the break- 
down of the semiclassical picture of electrons following 
classical trajectories that is valid at lower B. By contrast, 
no breakdown of the semiclassical picture has been found 
in previous work on large arrays of antidots Q because 
in those devices |j] the constrictions were much wider, 
each transmitting many modes. Our quantum simula- 
tions also reproduce the experimentally observed abrupt 
appearance of additional AB resonances as the antidot 
geometry is varied and the subsequent re-emergence of 
the regular sequence of conductance peaks; an explana- 
tion of this phenomenon is also suggested. 

II. THE MODEL 

The model used in our quantum simulations is a tight- 
binding Hamiltonian on a square lattice 

(1) 



l 



Here m and n label sites in the x and y directions in 
Fig.[l], t = h 2 /(2m*d 2 ), a = eBd 2 /Ti (the magnetic phase 
factor e ma is associated with hopping in the x direction 
since we use the Landau gauge A = (—By, 0,0) ||), B 
points in the z direction, d is the lattice parameter, — e 
is the electron charge, m* is the electron effective mass. 
We ignore the small Zeeman splitting between spin up 
and down. The electron potential energy (the site en- 
ergy W mn in equation ([!])) was assumed to be parabolic 
near the edges of the gates that define the boundaries of 
the wire and the antidots, and flat elsewhere. W mn was 
modeled as a sum of partial electron potential energy 
functions W g (u) associated with the individual gates g, 
with W g (u) = E F (u - a(l + s g )) 2 /a 2 for u < a(l + s g ) 
and W g (u) = for u > a(l + s g ). u is the distance 
from the edge of the gate, a is a scale of length that was 
chosen to be 0.05/i. s g is the dimensionless width of the 
depleted region around the gate. (In this work s g = 1 
for both gates that define the edges of the quantum wire 
and s g = Sd for both gates that define the antidots). 
Ep = Txnh 2 /m* is the electron Fermi energy measured 
from the bottom of the band and n = 3.47 x 10 -15 m -2 
is the 2D electron density far from the gates. As in the 
experiments, the distance between the outer gates that 
define the quantum wire in Fig.l was lfi, each antidot 
had a lithographic diameter of 0.2/i, and the spaces be- 
tween the two antidots and between the antidots and 
outer wire gates were 0.2/j,. In our simulations 100 lattice 
sites spanned the 1/i width of the wire. For a typical value 
of B = 0.25T this implies 5 lattice spacings d per mag- 
netic length {Ti/ (Be)) 1 / 2 . Thus although our numerical 
results reproduced the experimentally observed behavior 
of the device quite well, it should be emphasized that the 
model of the system used was a moderately coarse tight 
binding lattice rather than a continuum B . 



III. QUANTUM SIMULATIONS AND 
COMPARISON WITH EXPERIMENTS 



The electron wave functions and the multichannel 
quantum transmission and reflection matrices for the 
above model were calculated numerically using a stabi- 
lized transfer matrix technique [fiof . The two-terminal 
conductance G was then obtained from the Landauer for- 
mula 0] G= ^Tr(ttt). 

The calculated conductance G of the device is shown 
in Fig for two values of the width of the depleted re- 
gions around the antidots. Plot (b) [(c)] is for Sd = 2.050 
[1.850] for which the individual conductances G\ of the 
two constrictions between the antidots and the outer 
gates are both accurately quantized to 2e 2 /h [4e 2 //i]. 
I.e., one [two] transverse modes are transmitted perfectly 
through each constriction. The constriction between the 
antidots is pinched off. At the upper and lower ends 
of the range of B shown, the conductances of the two 
parallel constrictions add (as at B = |E2|) and the 



conductance of the device is close to G — 4e 2 /h [8e 2 /h]. 
However, in between G is depressed; this is the dip in the 
conductance of the device that is observed experimentally 
H due to the trapping of electrons in the orbit around 
the antidots. This effect is remarkably strong - near the 
center of the dip the calculated conductance falls to little 
more than half of its maximum value. As in the exper- 
imental data, the calculated G exhibits AB oscillations 
superimposed on the conductance dip, and as in the data 
the pattern of oscillations is much simpler when only one 
mode is transmitted through each constriction. (The ex- 
perimental trace (a) is in the single mode regime as is 
plot (b)). The conductance dip in the simulations is at a 
somewhat higher value of B than in the experiment; we 
did not tune the model parameters for a perfect match. 
This should be remembered when comparing simulations 
with experiment in Fig|| and in the top panel of Figjl]. 

A feature of our numerical results is that the AB res- 
onances (when they are strong) usually take the form 
of sharp peaks and smooth minima in the conductance; 
in the simpler case (b) the peak conductance equals the 
ideal value 4e 2 //i for perfect transmission through both 
constrictions. This means the AB resonant states give 
rise to resonant transmission. Trace (a) is consistent with 
what would be obtained by broadening a series of peaks 
with sharp maxima and smooth minima like those in plot 
(b); no kT broadening is included in plots (b) and (c). 

The locations in B of the peaks in the calculated con- 
ductance G are shown in Fig.|(a) for various widths Sd of 
the antidot depletion regions. The results in Fig.||(a) are 
in good qualitative agreement with the corresponding ex- 
perimental data plotted against the antidot gate voltage 
V g in Fig.|^(b). In both the experiment and simulations 
the pattern of conductance maxima is simple and regular 
when one mode is transmitted through each of the two 
constrictions (G\ = 2e 2 /h). But in the regime where 
two modes are perfectly or partly transmitted there are 
anomalies - a conductance maximum disappears as Sd in- 
creases near (sd — 1.8, B = 0.216T) and as Sd decreases 
near (sd = 1.97, B — 0.30T). But as Sd is varied further 
the remaining conductance peaks form into a periodic 
pattern once again. I.e., these anomalies resemble dislo- 
cations in a crystal. The same two dislocations are seen 
at somewhat lower B in the experimental data when two 
modes are transmitted through each constriction Jl3| . A 
series of experimental conductance traces around the dis- 
location at B = 0.237T is shown in the inset of Fig.||(b). 
Note that in the dislocation found experimentally on the 
high [low] B side of the trapping dip in the conductance, 
the conductance peak disappears as the antidot gate volt- 
age becomes more [less] negative. This is in agreement 
with the orientation of the corresponding dislocations in 
Fig.^(a) since increasing the depletion width Sd in the 
simulations means making V g more negative. 

Since the dislocations in the pattern of AB resonances 
are found in our simulations that are based on a model 
of non-interacting electrons, it seems unlikely that they 
are due to electron-electron interactions. However, both 
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experimentally and in the simulations the dislocations 
occur when more than one mode is transmitted through 
the constrictions. Indeed more complex dislocation pat- 
terns are seen in FigJ|(a) for Sd < 1.7 where a third 
mode is transmitted. More than one mode penetrating 
the constrictions should result in more than one trapped 
electron orbit of the antidot pair being present. This 
is consistent with the beat-like modulation of the am- 
plitude of the AB resonances visible in plot (c) of Fig.[l] 
(two modes transmitted), but absent in plot (b) (one 
mode penetrates). Our numerical study of the trans- 
mission probabilities associated with individual scatter- 
ing channels that contribute to the total conductance of 
the device showed that when two modes are transmitted 
through the constrictions the resonant structures in the 
different channels fall into two groups. The peaks in the 
transmission belonging to these two groups move out of 
phase with each other in the vicinity of the dislocations 
and cancel there resulting in the disappearance of a con- 
ductance peak. Thus it appears that destructive interfer- 
ence of out of phase resonances associated with at least 
two different orbits is responsible for the dislocations. 

The spacing of successive conductance peaks in B (the 
AB period AB) found in the quantum simulations is 
shown by the large open [full] circles at the top of Fig.[l] 
for the same Sd as in conductance plot (b) [(c)]. AB 
decreases with increasing B at low B but is station- 
ary at the conductance dip. (The anomaly in AB at 
B = 2.8T for case (c) is due to a dislocation whose core 
is at (s d = 1.65, B — 3.0T) in Fig.|(a)). Since AB could 
only be measured in a relatively narrow range of B at 
the conductance dip, the fact that the simulations show 
AB to be stationary there accounts for the independence 
of AB of B that was observed by Gould et al. The 
small symbols in the top panel of Figjl] show the mea- 
sured AB for several values of V g . The data shows some 
scatter but there is overall agreement with the quantum 
simulations - the measured AB also decreases with in- 
creasing B at low B and levels out at the conductance 
dip. Notice that overall AB is insensitive to Sd in both 
Fig.|l| and Fig.||(a), consistent with the insensitivity to 
the antidot gate voltage found experimentally 0. 

IV. ANALYSIS: QUANTUM SIMULATIONS VS. 
CLASSICAL SIMULATIONS AND FEYNMAN 
PATH INTEGRAL THEORY 

The fact that AB decreases with increasing B at low 
B may at first seem surprising since the area enclosed by 
a cyclotron orbit in free space decreases with increasing 
B. However, at low B the diameter of a cyclotron or- 
bit is larger than the distance between the gates defining 
the edges of the wire, so that a closed circular orbit en- 
closing the antidots can not occur. But the constrictions 
between the the gates and the antidots are asymmetric; 
one side is curved, the other straight - See left inset, 
Fig.ll]. This results in the electron wave, as it emerges 



from the constriction, being collimated at an angle of 
about 7r/6 relative to the edge of the straight gate [Q so 
that a closed non-circular quantum orbit of the antidots 
can form. 

A similar collimation effect also occurs classically [ pT[ , 
making closed classical low B orbits possible. We have 
performed classical calculations of such closed orbits for 
the same model electron potential as was used in our 
quantum simulations and calculated the AB phase as- 
sociated with these orbits using the semiclassical result 
h$AB = § m * v -^l ~~ eA e B. Here A e is the area enclosed 
by the orbit. The semiclassical AB period thus obtained 
for a closed classical orbit is shown by the line in the top 
panel of Fig.Q; the semiclassical AB period decreases with 
increasing B as in the quantum simulations but does not 
become stationary ]lq] . Since the same electron potential 
energy function was used in both the quantum and clas- 
sical simulations, the fact that the stationary AB period 
was only found in the quantum simulations is significant: 
It implies that the constant AB period observed experi- 
mentally is an inherently quantum phenomenon that sig- 
nals the breakdown of the semiclassical picture of elec- 
trons following classical trajectories that is valid at lower 
B. 

One can understand the effect of the confinement of 
the electron orbit by the edges of the quantum wire on 
the AB period at low B by considering a simplified model 
of the device and a class of Feynman paths depicted in 
the upper right inset of Fig.[l]. The edges of the wire 
are represented by infinite hard walls and the antidots 
by a narrow infinite barrier. In the symmetric gauge 
A = (—By/2,Bx/2,0), we parameterize the Feynman 
path by x = k{v/u>) cos(wt), y = (v/u>) sin(wt), where 
t is the time, v is the Fermi velocity and u> — v/r. r 
is half of the distance between the constrictions, k is 
the parameter that defines the shape of the particular 
Feynman path. The action for a closed path is then 
given by S k = § Ldt = Tm*v 2 (k 2 - 2kuj c /u + l)/4 
where T = 2ir/uj and lu c — e\B\/m*. The action is 
stationary with respect to the shape parameter k for 
k = k = lu c /lu. Thus the stationary Feynman path is 
circular as expected when lu c = u>, i.e., when the cy- 
clotron radius equals r. For lower B , oj c < uj and there- 
fore k < 1 and the stationary Feynman path is an ellipse 
whose area increases as the B increases. This suggests 
that the AB period should decrease with increasing B, 
qualitatively in agreement with our low B simulation re- 
sults. On closer inspection, this turns out to be correct 
- The propagator U(r, t, r, t + T) for the set of Feynman 
paths enclosing the barrier is dominated by e lS ' k ^ h , and 
therefore the AB phase is given (up to a constant) by 
4>ab « SJfi = Tm*v 2 {\ - uj 2 /lo 2 )/(AK). The AB pe- 
riod is then given by 2n — \A4>ab\ ~ AB \d4>AB / dB\ or 
AB w Ahm*TTuj 2 /(BTv 2 e 2 ). I.e., AB decreases as B in- 
creases. More general path integral calculations that con- 
firm this qualitative result will be presented elsewhere. 

At higher B the cyclotron diameter becomes smaller 
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than the width of the quantum wire, so that the confine- 
ment of the AB orbit by edges of the quantum wire should 
be less important. Thus at higher B one should expect 
the area of the quantum AB orbit to shrink and the AB 
period to increase with increasing B. The crossover to 
this higher B regime should occur when the cyclotron 
diameter is equal to the distance between the two con- 
strictions at the edges of the device, i.e., within the dip 
in the conductance of the device that is due to the trap- 
ping of electrons in the orbit encircling the antidots. This 
explains why the AB period in our quantum simulations 
begins to increase with increasing B at the high B edge 
of the resistance peak, as can be seen in Fig.|l|. Thus 
AB was effectively independent of B in the experimen- 
tal data of Gould et al. [|| because that data was taken 
at the crossover from the low B regime to the high B 
regime, which coincides with the trapping dip in conduc- 
tance 

As a further test of this explanation we carried out 
quantum simulations of transport in a similar device in 
which the straight edges of the quantum wire are further 
from the constrictions and should have less effect on the 
AB orbit at the lower values of B. This device is depicted 
in the lower right inset of Figj^. As expected from the 
above argument, in our simulations of this device the AB 
period increased with increasing B throughout the range 
of B in which AB conductance oscillations were present. 
Experiments by Hirayama and Saku jl7) on such a device 
also found the AB period to increase with increasing B, 
lending further support to the above theory. 



V. CONCLUSIONS 

In conclusion, we have shown that the observed 
Aharonov-Bohm resonance patterns of an antidot clus- 
ter can be accounted for by quantum simulations and 
have demonstrated that the previously unexplained con- 
stant period of the resonances signals the crossover from 
a classical ballistic regime at low magnetic fields to an in- 
herently quantum regime at higher fields. The crossover 
from classical to quantum ballistics should be an inter- 
esting new field of study. 
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FIG. 1. Two-terminal conductance G and AB period AB vs. B for device in left inset. Plot (a) [(b)]: experiment [quantum 
simulation] with one mode transmitted through each constriction. Plot (c): quantum simulation with two modes transmitted. 
Large open [full] circles: AB for case (b) [(c)]. Small open [full] symbols: experimental AB in one [two] mode regime. Solid line 
in top panel: semiclassical AB vs. B for closed classical orbit of antidot pair; s,i — 2.0. Left inset: Schematic of quantum wire 
with two antidots. Black and shaded regions are depleted of electrons. Upper right inset: Simplified model used in analytic 
path integral theory. Lower right inset: Device in which the straight edges of the wire are further from the constrictions. 

FIG. 2. Positions in B of conductance maxima obtained from (a) quantum simulations vs. width Sd of the antidot depletion 
regions and (b) measurements vs. antidot gate voltage V g . Quantized plateaus in conductance Gi of individual constrictions 
are indicated. Curves are guides to the eye. Inset: Measured conductance G [for various V g from -1.16V (top trace) to -1.32V] 
vs. B (from 0.21T to 0.26T) around dislocation at V g = -1.22V, B = 0.237T. Traces offset for clarity. 
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